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We study the dynamics of the travelling interface arising 
from a bistable piece-wise linear one-way coupled map lattice. 
We show how the dynamics of the interfacial sites, separating 
the two superstable phases of the local map, is finite dimen- 
sional and equivalent to a toral map. The velocity of the 
travelling interface corresponds to the rotation vector of the 
toral map. As a consequence, a rational velocity of the trav- 
elling interface is subject to mode-locking with respect to the 
system parameters. We analytically compute the Arnold's 
tongues where particular spatio-temporal periodic orbits ex- 
ist. The boundaries of the mode- locked regions correspond 
to border-collision bifurcations of the toral map. By varying 
the system parameters it is possible to increase the number 
of interfacial sites corresponding to a border-collision bifurca- 
tion of the interfacial attracting cycle. We finally give some 
generalizations towards smooth coupled map lattices whose 
interface dynamics is typically infinite dimensional. 

I. INTRODUCTION 

Coupled map lattices (CML) where introduced by 
Kaneko [1983, 1984] as a paradigm for the study of 
spatio-temporal complexity such as turbulence [Kaneko, 
1986, 1989, Beck 1994], convection [Yanagita & Kaneko, 
1993], open flows [Willeboordse & Kaneko 1995, Kaneko 
1985], patch population dynamics [Hassell et al, 1995, 
Sole & Bascompte 1995], etc. An interesting spatio- 
temporal feature is the emergence of coherent travelling 
structures [Kaneko 1992, 1993] along the lattice. In this 
paper we use a simple example of a piece- wise linear one- 
way CML with two superstable basins in order to restrict 
the front dynamics to a finite dimensional system. At the 
beginning of Sec. 2 we recall some features in a particular 
region of the parameter space where the travelling inter- 
face consist of a single site. The dynamics in this case 
can be reduced to a one-dimensional map of the circle. 
The velocity of the travelling interface is computed by 
means of the rotation number of the circle map and thus 
the mode-locking regions of the velocity with respect to 
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the system parameter are computed analytically. Later 
in Sec. 2 we perform an equivalent approach for the case 
of a 2 sites travelling interface. Here the dynamics of the 
travelling interface can be reduced to a two-dimensional 
toral map whose rotation vector gives the travelling ve- 
locity. We give a method for computing the mode-locking 
regions in this case. In Sec. 3 we present the Arnold's 
tongues for higher dimensional interfaces and show the 
bifurcation diagram as the number of sites in the inter- 
face is increased. Finally we give some comments on the 
generality of our results for generic bistable local maps 
whose attractors do not possess superstable basins and 
thus their corresponding travelling interface is typically 
infinite dimensional. 

A CML is a dynamical system with discrete space, dis- 
crete time and continuous state space. We could think 
of CMLs as a generalization of cellular automata whose 
state space is discrete [Chate & Manneville, 1989, 1990]. 
Consider a one-dimensional collection of cells labeled by 
the integer index i and with dynamical variable Xt(i) at 
time t. The lattice could be infinite or finite with peri- 
odic or fixed boundary conditions. For the present study 
the particular choice of lattice is not crucial since we deal 
with the propagation of finite dimensional localized struc- 
tures with spatio-temporal periodicity. Therefore, for the 
numerical experiments, it suffices that the lattice is large 
enough to enclose the localized front for a whole period. 
The CML dynamics consists of two independent stages: 
the local dynamics and the coupling dynamics. The for- 
mer is the application of the one-dimensional local map 
fi(x) to every site and the latter couples the dynamics 
by means of a weighted sum over a neighbourhood i>i\ 

x t ii) local fi(xt(i)) 

> (1) 

coupling x t+1 (i) = e v J Vi (x t {i)), 

where the coefficients e Vi determine the type of coupling 
interaction. A physically meaningful interaction will have 
typically a limited range, with decreasing strength for 
distant neighbours. The local map and neighbourhood 
are often taken to be the same for all sites (/; = / 
and vi = v) indicating an homogeneous dynamics - 
homogeneous CML. Also one asks that ^2e Vi = 1 as 
a conservation law, since failure to do so may lead to 
non-boundedness of the state as time tends to infinity. 
The two most widespread models of CML are 
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x t+1 (i) = (1 - e)f{x t {i)) + - (f(x t (i - 1)) + f(x t (i + 1))) 

(2) 



and 



x t+1 (i) = (1 - e)f{x t (i)) + ef(x t (i - 1)), 



(3) 



which are called diffusive and one-way CML respectively. 
The coupling parameter e is set to satisfy < e < 1. 

In this paper we study interfaces separating two dif- 
ferent phases in bistable one-way CMLs. For two differ- 
ent phases to coexist in the lattice we consider a local 
map with two and only two stable fixed points x*_ and 
x* + (x*_ < x+), i.e. a bistable local map, then each of 
the two homogeneous states Xt(i) = x± is stable for the 
CML dynamics. We now consider minimal mass states 
of the lattice whose sites are arranged in non-decreasing 
order (x t (i + 1) > x t (i)), and choose local maps that 
are continuous and non-decreasing. For such maps it is 
possible to have an interface, sites separating the stable 
x± phases, that travels along the lattice [Carretero et 
al. 1997a, 1997b]. Travelling fronts in generic bistable 
CMLs are typically infinite dimensional since the attrac- 
tion towards the stable points x± is an infinite process. 
Nevertheless the convergence towards x± is exponential 
(for stable fixed points). This exponential convergence 
ensures a localization of the front despite the fact that 
the interface involves an infinite number of sites. 

The simplest of all the local maps with two superstable 
regimes is the family of piece-wise linear local maps f a 
defined by 

— 1 if x < —a 
fa(x) = { I x if -a < x < a < a < 1, (4) 
1 if x > a 

and by the step function between —1 and 1 centered 
at the origin when a — 0. The fixed points x*_ = —1 
and x* + = 1 are superstable, with superstable basins 
5- = [— 1,— a] and S+ = [a, 1] respectively, and thus 
the homogeneous phases x(i) = x± are superstable as 
well. The parameter space is now (e,a) 6 [0,1] x [0,1]. 
The piece-wise linear family f a possesses the important 
property of collapsing orbits that get (1 — a)-close to x±. 
The propagating front is then subject to a cut-off near 
Xj- ensuring a finite number of interfacial sites (falling 
in the linear regime [—a, a], for < a < 1). Therefore, 
instead of dealing with an infinite dimensional interface 
— that is the case for generic CMLs — we are let with 
a finite number of interfacial sites allowing us a precise 
description of the front dynamics. 



II. ONE- AND TWO-DIMENSIONAL 
INTERFACE DYNAMICS 



A. One-dimensional interface dynamics 

In a previous paper [Carretero et al., 1997a] we showed 
how the dynamics of the whole lattice could be reduced to 
a one-dimensional map for a particular region of the pa- 
rameter space. We now briefly recall some of the results 
in order to use the same approach for the iV-dimensional 
case. 

Consider a general minimal mass state with N sites in 
the interface. The interface, separating the two phases, 
is well-defined in the case of the piece-wise linear map 
f a since any site falling into the superstable regions S± 
is mapped to ±1 by f a in one iteration. Therefore we 
define the interface as the collection of sites belonging 
to the unstable region U = [—a, a]. We will often use 
the term interface for denoting the interval U itself. The 
general form of a minimal mass state is 

X t = (..., -l,-l,x t (k + l),...,x t (k + N), 1,1,...), 

with |xt(i)| < a for i = k + 1, . . . , k + N. Without loss 
of generality one may choose k — 0. The image of X t by 
the one-way CML (|^) is 



X t +i = (...,-l,g(-a,xi),g(xi,x 2 ),. 

g(x N -i,x N ),g(x N ,a), 1, . . .), 



(5) 



where now the explicit time dependence is dropped and 
the subscript denotes the label of the site. The function 
g, defined by 



. . 1 - e e 

g{xi,x 2 ) = x 2 + - xi, 

a a 



(0) 



originates from the combination of the local dynamics 
and the coupling. At both extremes of the interface the 
function g only uses explicitly one site since the other is 
— 1 or 1, left and right extremes respectively. By defining 
now the maps 



fo(x) = g{-a,x) 



1-e 



x — e 



fi(x) = g{x,a) = | x + (1-e), 

it is possible to write (Eh as 

X t +i = (...,~l,f Q (x 1 ),g(x 1 ,x 2 ),... 

g(x N -i,x N ), f!(x N ),l, . . .), 



(7) 



(8) 



where now the map /o (/i) originates from the combi- 
nation of the local map and the interaction between the 
left-most (right-most) site of the interface with the neigh- 
bouring homogeneous region, while the map g originates 
from the interaction between interfacial sites. 
By defining 

a(l — a — 2e(l — ej) 

7(£, a) = 7 - - 7+ = 7, ^ (9) 

e(l -e) 
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where 

a(e — a) a(a + e — 1) 



it is possible to show [Carretero et al., 1997a] that if 
7 > the state of the whole lattice is reduced to the 
one-dimensional circle map 





' h(x) 


if 


x G [—a, 7 




a 


if 






fc fo(x) 


if 


x G [7- , a] 



where T — (7+, 7-) and its length (with sign) is given 
by (^J) and is referred as the gap. The 2-parameter map 
of the circle $ e , a , called auxiliary map, accounts for the 
dynamics of the only interfacial site in the non-negative 
gap case. The auxiliary map maps the interfacial site 
throughout the evolution of the interface and the posi- 
tion of the latter is given by a symbolic dynamics repre- 
sentation. 




—a a —a a 



FIG. 1. Moments of birth a) and death b) of the v = 1/3 
tongue for a — 0.5. The v = 1/3 tongue begins (ends) when 
the interfacial orbit, mapped through the auxiliary map <& £ , a , 
collides with the left-most (right-most) extreme of the gap. 

The auxiliary map is a map of the circle and it is pos- 
sible to show [Carretero et at, 1997a] that the velocity 
of the travelling interface is given by its rotation num- 
ber. The superstable region (the gap) of the auxiliary 
map induces a mode-locking of the rotation number as 
the parameters are varied. All the mode-locking regions, 
Arnold's tongues, where any particular rational velocity 
< v(e) = p/q < 1 exists, may be computed analyti- 
cally by finding the (e, a)-interval for which the gap in- 
terval intersects the identity line after q iterations passing 
p times through the upper region of the auxiliary map. 
The extremes cases for this to happen arise from a border- 
collision bifurcation [Maistrenko et al., 1995] when an or- 
bit starting in the gap falls at the beginning (end) of the 
gap interval after q iterations. As an example we show 
in Fig. [I] the v = 1/3 plateau limit orbits for a — 0.5. 



B. Two-dimensional interface dynamics 

When the gap size is non-negative the whole lattice 
was reduced to the one-dimensional auxiliary circle map. 
For negative gap size there are more than one site in the 
interface at the same time. Thus a one-dimensional re- 
duction of the whole lattice does not seem attainable. 
However, thanks to the superstable basins S±, the num- 
ber of sites in the interface is still finite — less than or 
equal to N — and one could reduce the dynamics of the 
whole lattice to a A-dimensional system containing the 
information of the interface sites. 

1. Minimal N -state layers 

Before reducing the dynamics, one has to find the pa- 
rameter regions where a particular number of sites is 
present in the interface. Define a minimal mass N -state 
or minimal N -state to be a minimal mass state having, 
during its evolution, a minimum of A — 1 sites in the 
interface. The spatial discretization induces a change in 
the number of sites in the interface, so that the number 
of sites in the interface for a minimal TV-state cannot be 
constant, it varies during its evolution between A — 1 and 
A. 

Suppose that the front shape at any time t is given by 
xt(i) = h(i — vt), where v is its travelling velocity and h is 
a real valued function independent of t. The sites falling 
in the interface, are the sites contained in the interval 
H = (/i _1 (— a), h (+a)). Since the length of H and the 
distance between nodes remains constant, the number of 
sites in H can only vary by one during the evolution. In 
the case of a minimal A-state the number of sites in the 
interface is then A or A— 1. Nevertheless it is possible to 
remove this duality by a simple trick as it is shown in the 
next section. Using the same notation, the non-negative 
gap case corresponds to a minimal 1-state. 

In order to show that a minimal A-state possesses A 
or A — 1 sites in the interface we supposed that the front 
shape is given by a function h. If the interfacial orbit 
is periodic it is easy to construct the front shape h by 
considering any curve h(i — vt) passing through all the 
points of the orbit. It is interesting to notice that such 
h is not unique, in fact there exists an infinite number of 
possible /i-functions passing through the finite collection 
of required points. On the other hand, when the velocity 
is irrational the interfacial orbit consists of an infinite 
number of distinct points. It is possible to reconstruct the 
front shape for irrational v by superimposing snapshots 
of the front shape taken in a co-moving reference frame 
with velocity v. Extensive numerical experimentation 
suggests that the travelling interface is indeed given by an 
increasing function. Further reconstructions of travelling 
interfaces in more general types of CMLs are consistent 
with this result [Carretero et al., 1997c]. 
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FIG. 2. The v = 1/2 tongues in the minimal 2-state layer. 
The tongues a) and b) are obtained by solving the inequalities 
(pq) corresponding to the two possible period-2 orbits. 
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FIG. 3. Tongues in the TV- dimensional minimal TV-state 
layers. The tongues for the velocities v = 1/2, 2/5, 1/3, 2/7 
and 1/4 (from right to left) in the minimal iV-state layers for 
N = 1, . . . , 8 are presented. 

In Fig. H we depict the (e, a) region where the minimal 
2-state is present. The boundary between the minimal 
1-state and the minimal 2-state layers is the zero-gap 
curve. Above that curve, the gap is negative, and there 
appears to be an infinite family of layers, labeled by JV, 
and separated by boundary curves (cf. Fig. ||). 

2. Two-dimensional auxiliary map 

Throughout this section consider the case N = 2, 
i.e. take (e, a)-parameter values lying in the minimal 2- 
state layer. One would like to derive an auxiliary map 
that contains all the information of the interface dynam- 
ics as the auxiliary map $ 6iQ does for the non-negative 
gap case. In order to find such map one has to take into 



account all possible evolution combinations of a minimal 
2-state. The general form for a minimal 2-state X t is 

X t = (-1, xi, x 2 , 1, ■ • ■), 

but because the size of the interface may be 1 or 2 one 
has the two possibilities —a < x\ < x 2 < a or —a < x\ < 
a < x 2 < 1. Therefore, the first site, x±, is always in the 
interface while the second, x 2 , may or may not be in the 
interface. If the second site is not in the interface one 
could assign to it any value between +a and +1 and its 
image by the piece-wise linear map f a would remain the 
same. In this case we choose to reduce a < x 2 < 1 to x 2 — 
a and the dynamics is not altered. Using this reduction 
the site x 2 is now always included in the interface and it 
would be possible to remove the duality of having 2 or 1 
sites in the minimal 2-state layer. 




FIG. 4. Example of the two-dimensional auxiliary 
map Q,(xi,X2)=(ftl(xi,X2),fl2(xi,X!i)) accounting for the 
interfacial sites dynamics in the minimal 2-state for 
(e,o) = (0.365,0.7). 

It is possible to prove that the pair of interfacial sites 
(xi, x 2 ) is mapped to (a;^, x' 2 ) in 6 different ways depend- 
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ing on their value. The dynamics then induces a two- 
dimensional map 0(0:1,^2) = (%'i,x 2 )- An example of 
the map f2, obtained by considering all the possibilities 
is depicted in Fig. || by separating it into its two compo- 
nents x[ — Q,i(xi,X2) and x' 2 — 0,2(2 1, X2). The shaded 
areas in Fig. |] correspond to unreachable situations for 
a minimal 2-state. Since X t is a minimal mass state, X2 
has to be larger or equal to x\, thus we eliminate the area 
where X2 < x\ — light shaded area. On the other hand, it 
is straightforward to prove that if x\ and X2 are both at 
the same time in the interval [7-, 7+], the succesive state 
is a minimal 3-state and thus we eliminate this possibil- 
ity — dark shaded area. These forbidden areas are never 
reached by the dynamics. 

For the 2-dimensional case (minimal 2-state layer) the 
auxiliary map £1 is a toral map in two dimensions — a 
toral map is the generalization of a circle map to more 
than one dimension (cf. [Baesens et ai, 1991]). Gener- 
alizing the idea of rotation number in more dimensions 
one could think of a rotation vector [Baesens et al., 1991] 
whose entries correspond to the rotation number in every 
component. The rotation vector p = (p\, p2) of f2(ii, X2) 
is then two-dimensional. But since both sites (x±, X2) be- 
long to the same interface, that is moving with a defined 
velocity, the two components of the rotation vector have 
to be equal (p\ = pi)- In other words, the site X\ has to 
travel at the same velocity than the site X2 in order the 
interface to remain in the minimal 2-state layer. The only 
way X\ would have different rotation number than X2 is 
for an interface that is increasing (or decreasing) in size. 
Therefore instead of taking the whole rotation vector one 
may use a single scalar to describe the rotation around 
Q,(xx,X2)- Thus, the velocity of the travelling interface 
in the minimal 2-state layer is given by this scalar, which 
will be simply called from now, the rotation number of 
^(xi, re- 



values of (e, a) [Carretero et ai, 1997a]. However, in the 
two-dimensional case, there are different possible evolu- 
tions for a given velocity. This is due to the fact that 
with more sites in the interface a wider selection of or- 
bits is possible. The two orbits of a minimal 2-state giv- 
ing a velocity v = 1/2 correspond to a) X2 > a all the 
times and b) alternating one and two sites in the inter- 
face. First of all the orbit has to be periodic (period two). 
Therefore we must have a) fi(x' 1 — fo(x\)) = x\ and b) 
g(x'i = ft)(xi),x' 2 = fi(xi)) = xi, i.e., 



b) 



x x 



a(a(l 
•2 



e)-e 2 ) 



-e(l-e) 

n(l - 2e) 
J + 2e(e-l)' 



(11) 



Next, one has to verify that the sites fall in the right 
intervals. That is 

a ) /l(^l) > a and fo(%'i = fo(xi)) < -a, 

b ) fo(x'i = fo(xi)) < -a and fi(x' 2 =f 1 (x 1 ))>a. 

(12) 



Combining Eqs. (11) and ( pL2|) gives the conditions that 
e and a must satisfy in order to have the period-2 orbits. 
In Fig. ^ we display the regions where these conditions 
are satisfied in the minimal 2-state layer. 

The same procedure may be applied to find the tongues 
for any rational velocity in the minimal 2-state layer. For 
every chosen < v = p/q < 1 there are q possible dif- 
ferent orbits, by combining 1 and 2 sites in the interface, 
and each one has a corresponding tongue in the minimal 
2-state layer. 



3. The tongues for the two-dimensional case 

The two-dimensional auxiliary map f2 is piece- wise lin- 
ear — it is a combination of planes — and possesses the 
analogue of the gap T for the one-dimensional case: the 
region [— a, 7_] x [7+, a], see dashed plateau in Fig. ||. 
This region, which we refer to as the two-dimensional 
gap and denote by T2, acts in the same way as its one- 
dimensional analogue. Any (x\, rz^-orbit falling into T2 
is superstable and therefore is parametrically stable to 
perturbations. This parametric stability gives rise to the 
mode-locking of the velocity. Thus one expects the mini- 
mal 2-state layer to contain Arnold's tongues in the same 
way the minimal 1-state layer does. 

Let us derive the tongues for v = 1/2 in the two- 
dimensional case. First of all one has to choose the orbit. 
In the one-dimensional case the orbit was uniquely de- 
termined by a symbolic coding of the velocity for fixed 



III. iV-DIMENSIONAL INTERFACE DYNAMICS 

Let us now consider the case of a minimal iV-state layer 
for N > 2. The dynamics in that layer may be reduced 
to a A^-dimensional auxiliary toral map that controls the 
iV-tuplet of sites in the interface. Again, thanks to the 
superstable region of the local map, this A^-dimcnsional 
toral map will have a A-dimensional gap Tn such that 
any Af-tuplet in maps one of its components to a. The 
parametric stability is again induced by the presence of 
the gap r^v when the latter exists. 
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FIG. 5. Successive states of the lattice in the u = 1/3 
tongue. Every time a site touches the boundary of the su- 
perstable region, in order to pass to a new configuration, the 
tongue e-width is zero. The states from a) to i) correspond to 
increasing values of a inside the v = 1/3 tongues. The state 
a) is a minimal 1-state while the state b) is the transitional 
state to a minimal 2-state (zero gap). The states c), e) and 
g) correspond to the three different possibilities for a mini- 
mal 2-state and d) and f) correspond to the transition points 
between these three states. State i) is already is a minimal 
3-state and h) is the transition point between the minimal 
2-state to the minimal 3-state. 

Thus, in every minimal iV-state layer there are q 
tongues corresponding to the velocity v = p/q given by 
the q possible combinations for the periodic-g orbit of 
the iV-dimensional auxiliary map. The regions where 
each tongue exists are given by a system of inequalities 
that the iV-tuplet has to satisfy for the orbit to undergo 
the right combination of N and N — 1 sites in the in- 
terface. In Fig. H we show the tongues for the veloci- 
ties v = 1/2, 2/5, 1/3, 2/7 and 1/4, computed numer- 
ically, in the minimal JV-state layers for N = 1, ... ,8. 
We call the tongues in the minimal iV-state layer the 
sub- N -tongues since they emanate from the principal 
mode-locking tongues, the principal tongues, for the one- 
dimensional case. Therefore, from every principal tongue 
there are q corresponding sub-iV-tongues in each iV-layer. 
The structure of these sub-tongues is then self-similar 
and repeats itself in every layer as it may be observed in 
Fig. | 

It is interesting to notice that because of the conti- 
nuity of v(e,a) all the sub-A^-tongues touch each other. 
Furthermore, the sub-iV-tongues touch in a single point, 



i.e. in that particular point the width of the e-mode- 
locked plateau (for a fixed value of a) is zero. This phe- 
nomenon repeats itself all along each family of sub-JV- 
tongues and it happens every time a site of the interface 
touches the boundary of the superstable regions S± when 
varying the (e, a)-parameters in order to go from one of 
the q possible combinations of the interface orbit to the 
next. In Fig. ^ we show the state of the lattice at dif- 
ferent stages of the v = 1/3 tongue. The e-width of the 
v = 1/3 tongue is zero in the transitional case between 
two different interfacial orbits (cases b), d), f) and h)). 
One may consider this switch from one interfacial orbit to 
the next as a border-collision bifurcation of the attract- 
ing cycle of the-iV-dimensional auxiliary map. In order 
to illustrate when this bifurcation takes place we show in 
Fig. H the bifurcation diagram of the attracting cycle of 
the interface in the v = 1/3 tongue as the parameter a 
is varied. From the figure it is possible to observe that 
when a stable point of the interface touches the boundary 
of the superstable region (thick dashed diagonal lines a 
and —a) the attracting cycle changes because a further 
site is added to it. Where this happens (vertical dashed 
lines) the e-width of the tongue is zero (see Fig. ||). 
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FIG. 6. Bifurcation diagram of the attracting cycle of the 
interface sites in the v = 1/3 tongue family. We plot the 
stable periodic 3 orbit of the interfacial sites as a function 
of the parameter a after 10 000 transients. The values of e 
were taken such that the pair (e, a) remained in a path inside 
the v = 1/3 mode- locked tongue (see dashed line inside the 
tongues in Fig. |B|). 

IV. GENERALIZATIONS 

For every family of tongues for a particular mode- 
locked ratio there is an equivalent scenario as the one 
presented for v = 1/3. The mode- locking regions could 
be computed analytically in every layer along with the 
attracting cycles and bifurcation diagrams. This is made 
possible by the simple form of the piece-wise linear local 
map. In more general CMLs, one way or diffusive, with 
bistable smooth local maps, the mode-locked regions still 
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exist though they are not as large as in the piece-wise lin- 
ear case because they typically involve an infinite number 
of sites in the interface [Carretero, 1997b]. The reduction 
of the dynamics to a toral map for the piece-wise lin- 
ear local map requires a finite number of interfacial sites. 
Therefore it seems impossible to apply a toral map reduc- 
tion of the dynamics for a generic CML. Nevertheless, if 
there exists an invariant travelling interface h(i — vt) it is 
possible to reduce the dynamics of an infinite interface to 
a one-dimensional circle map whose rotation number de- 
scribes the velocity of the travelling interface [Carretero 
et al, 1997c]. 

The border-collision bifurcations corresponding to 
transitions between families of sub-tongues prevails for 
more general CMLs and a similar scenario as the one 
depicted here for non-negative gap still exists. On the 
other hand, the border-collision bifurcations correspon- 
ding to transitions between sub-tongues of the same fam- 
ily, i.e. with the same rotation number, do not occur for 
smooth local maps. For the piece-wise linear CML these 
bifurcations came from a collision of a site with the bor- 
der of the interface. For generic CMLs the interface is 
the whole interval [a;!H , and thus it is impossible to 
have this kind of border-collision bifurcations. 
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